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A general scenario for electronic charge migration in finite samples illuminated by an intense laser 
pulse is given. Microscopic calculations for neon clusters under strong short pulses as produced 
by X-ray free-electron laser sources confirm this scenario and point to the prominent role of field 
ionization by strong internal fields. The latter leads to the fast formation of a core-shell system with 
an almost static core of screened ions while the outer shell explodes. Substituting the shell ions with 
a different material such as helium as a sacrificial layer leads to a substantial improvement of the 
diffraction image for the embedded cluster thus reducing the consequences of radiation damage for 
coherent diffractive imaging. 
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The advent of sources producing short and intense pulses 
of light with frequencies ranging from soft to hard X- 
rays d.j opens a new parameter regime for light-matter 
interaction. Dynamics with a quick removal of many 
electrons from their bound states in the atoms of the 
sample irradiated has scarcely been explored before, ft 
is important to have a good insight into this dynamics 
in order to realize one of the goals of X-ray free elec- 
tron lasers (XFELs) which deliver such light pulses: the 
single-shot coherent diffractive imaging of finite samples 
with atomic resolution 0]. X-ray photons scatter elasti- 
cally from a single non-periodic molecule, like a protein, 
and are subsequently recorded on a CCD as a continu- 
ous diffraction pattern The estimated total photon 
flux required for single-shot diffraction imaging exceeds 
10 12 photons per (lOOnm) 2 and is at the limit 

of what XFELs can deliver. Above all, such high photon 
fluxes lead to radiation damage by multiple ionization 
and subsequent Coulomb explosion of the sample which 
degrades the diffraction image if not destroys it. 

For possible means to overcome this handicap a 
detailed understanding of the time-dependent multi- 
particle dynamics is required which is very interesting 
in its own right. Apart from simulations in the X-ray 
regime 0, H, 0, H, H, EH [H| , this problem has been ap- 
proached from the perspective of intense long- wavelength 
laser pulses. There, it was found that finite systems, such 
as atomic clusters or nano-droplets, display behavior fun- 
damentally different from single atoms or bulk systems 
[H, El • The large space charge of ions in a finite volume 
creates strong electric fields easily exceeding the strong 
external field of the laser. Similarly, very strong internal 
fields are generated by single-photon ionization in X-ray 
pulses since many photo electrons can leave the sample 
due to the high photon energy. In both cases the strong 
internal field leads to the release of additional electrons, 
mainly from the surface. In the X-ray regime this effect 
has profound consequences for diffractive imaging, as we 
will show. 

As long as the laser frequency is sufficiently high to 



trigger initial ionization, a unified scenario with four 
phases emerges for a sample illuminated by an intense 
light pulse: i) Electrons are photo-detached from the 
atoms by the laser and leave the cluster, ii) The sub- 
stantial ionic charge developed during the first phase 
generates a large electric field particularly at the cluster 
surface, leading to almost instantaneous field ionization 
of surface atoms. These electrons migrate towards the 
cluster center, iii) An electron plasma is established, ini- 
tially formed by field-ionized electrons. For sufficiently 
large samples secondary electrons (from Auger decay 
or electron-impact ionization) and even photo-detached 
electrons are trapped by the cluster potential, "feeding" 
the plasma, iv) The electron plasma shields the core of 
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FIG. 1: (color online.) Electron trapping and strong-field ef- 
fects of homogeneously charged neon clusters with an average 
charge q per ion and a cluster radius R. Trapping of elec- 
trons detached by 12 keV photons (green-hatched region) and 
Auger decays (blue-hatched region). The thick lower border 
(dark-green and dark-blue bands) take into account that elec- 
trons are released from any position in the cluster. Red lines: 
Ionization of surface ions at various charge states due to the 
internal radial field of the charged cluster estimated by the 
Bethe rule, see text. 
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the cluster so that only the (highly charged) outer shell 
explodes. 

This scenario occurs for many combinations of parame- 
ters (sample size and density, pulse length, frequency and 
intensity) provided that the quiver amplitude x = F/to 2 , 
i.e., the excursion of electrons induced by the oscillating 
laser field with amplitude F and frequency to, is much 
smaller than the characteristic size (linear extension X) 
of the sample, x <C X. In the X-ray regime @, Q this 
is fulfilled even for high field strengths, but it has been 
shown to occur also for much smaller frequencies in the 
VUV [HI ■ Recently, this electron migration has been el- 
egantly verified experimentally with an ad hoc desi gned 
core-shell system of a xenon-argon composite cluster [15] . 

Figure [1] shows at which combination of cluster radius 
and average charge per ion trapping of photo- and Auger 
electrons is expected for neon clusters exposed to a 12 
keV photon pulse, relevant for atomic-resolution imaging. 
One sees immediately, that large clusters trap even X- 
ray photo-electrons while samples of moderate size trap 
only the much slower Auger electrons. Field ionization, 
in contrast, occurs from the smallest sample sizes on. 
What contributes even more to their prominent role for 
the dynamical screening of ions is the fact that the field 
ionized electrons are not only available in a sample of ar- 
bitrary size but they are also the first electrons available 
for screening once the ionizing laser pulse sets in. Auger 
electrons, on the other hand, appear only about 2.5 fs 
after a K-shell photo-ionization event in neon, when the 
core hole undergoes an Auger decay. 

The critical field £ c for ionization of an atom can be 
estimated by the Bethe-rule £ c = E ip 2 /(AZ), where 
Ei p is the ionization potential and Z is the charge of the 
binding core. Hence, fields produced in neon clusters as 
small as R ps 10 A are sufficient for multiple ionization of 
the surface atoms (red lines in Fig.[T]) while larger clusters 
with radii R > 50 A reach the critical field strength at 
the surface with only a small fraction of photo-ionized 
atoms. This suggests field ionization as an extremely 
efficient and fast mechanism for an electron migration 
that leads to screening of the core of the sample. 

To assess the effect of field ionization quantitatively we 
have calculated the dynamics of neon clusters of N — 50 
(R ps 6 A) up to N = 15000 (R ps 42 A) atoms under in- 
tense, short (T = 1. . . 20 fs) pulses with a photon energy 
of E p h — 12keV focused to an area of A = (100 nm) 2 . 
The pulse amplitude in time is given by 



F(t) = /W' 



for \t\ < (ttT)/(2c) (1) 



with c = 2 arccos(2 -1 / 4 ) ps 1.144 so that T has the mean- 
ing of the full- width at half maximum pulse length of the 
intensity I(t) = F 2 (t). For a fixed number of n = 10 12 
photons in the pulse, the peak intensity follows as 
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FIG. 2: (color online.) Radial electric field (upper panels) 
and radial charge density (lower panels) for a Neisoo cluster 
and T = 10 fs pulse without the effect of field- ionization (left) 
and with field ionization (right). White lines show the average 
radial coordinate of cluster surface shell (the 5% outermost 
ions) and of cluster-core surface shell (the 5% outermost ions 
of the inner half of the cluster). 



Note that for different pulse lengths T, as discussed be- 
low, the peak intensity 7 pea k changes. Our approach is 
based on a classical molecular dynamics simulation sim- 
ilar to the one described in [H, [T3| , with photo-ionization 
and Auger decay treated by quantum-mechanical rates. 
The important effect of field-ionization is included by 
always propagating the least bound electron of each 
atom/ion as a classical particle in the field of all other 
ions and electrons. An ionization event is taken to have 
occurred, when this electron leaves its mother ion beyond 
a threshold radius. We have found the system dynamics 
to be insensitive to the particular choice of threshold ra- 
dius, so long as we stay within the natural limits of half 
the nearest neighbor distance from above and a minimum 
radius which avoids counting bound electrons with large 
excursion radius as ionized, despite their immediate re- 
turn to the allowed region. Within this range we have 
settled on the threshold radius value of 2.2 atomic units. 
Our approach has the advantage of incorporating the to- 
tal field of all charged particles, therefore including both 
electron-impact ionization and ionization due to static 
fields [lH. We made use of an implementation [l|| of 
the fast-multipole method [2(| to calculate the Coulomb 
interaction of electrons and ions in order to bring calcu- 
lation times for the larger systems to a manageable level. 

Paradigmatically, we first consider a Neisoo cluster un- 
der a 10 fs pulse. To clearly identify the effect of field 
ionization, we have performed two calculations for each 
parameter set. A full calculation including ionization 
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from quasi-static internal fields and a reduced calcula- 
tion without this effect. The latter amounts in our ap- 
proach to exclude the classical propagation of a bound 
2p electron with every atom/ion. At peak intensity, half- 
way through the pulse roughly each atom will have been 
singly photo-ionized. However, with a radius of R ss 20 A 
this cluster is not yet large enough to capture any of 
the photo-electrons, cf. Fig.[T] Without field-ionization 
this leads to a homogeneous charging of the cluster, see 
Fig. [2k The resulting electric field (Fig. [2k.) exhibits a 
linear increase inside the cluster and a 1/r dependence 
outside the cluster. The maximum field at the surface 
is about 1.5 atomic units which is sufficiently high for 
further ionization and should not be neglected. 

Field ionization generates many plasma electrons in 
the cluster. More specifically, at peak intensity there will 
already be an average of two plasma electrons per cluster 
ion created in this way. Thereby, the cluster effectively 
becomes a charged conducting sphere: The electrons neu- 
tralize the ionic charges at the center of the cluster and 
the excess positive charge is localized on the surface as 
can be seen in Fig.[5Ji in accordance with phase iii) de- 
scribed before. As anticipated, this highly efficient charge 
migration begins almost immediately with the photo- 
induced charging of the cluster and takes place on a sub- 
femtosecond time scale preceding the trapping of Auger 
electrons as K-shell holes decay much more slowly. The 
emergence of a neutral core with a positively charged 
shell leads to a fundamentally different electric field, as 
is shown in Fig.[5j3: Although the field at the surface is of 
similar strength as in the homogeneously charged cluster, 
there is now a nearly field-free region in the cluster center. 
Within this region the ion motion is therefore drastically 
suppressed at the expense of a violent explosion of the 
more strongly charged outer shell. 

The averaged radial trajectories of the surface ions 
starting at r « R and of the shell enclosing the inner half 
of the cluster ions starting at r a R/2 1 / 3, (white lines 
in Fig. [2]) indicate this twofold dynamics of the cluster 
ions. Clearly and in contrast to the reduced calculation, 
the cluster core in the full calculation displays almost no 
expansion throughout the pulse. 

The mean displacement 
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of the ions from their initial positions (t — > — oo) provides 
a systematic and more quantitative measure of the effect 
of charge migration and its dependence on pulse duration 
T and sample size N. It is shown in Fig. [3k. for various 
pulse lengths T (but a fixed photon number per pulse) 
in the full and the reduced calculation at t = 0, the time 
of peak intensity from which the largest contribution to 
the diffraction pattern in an imaging experiment can be 
expected. We find a strong amplification of ionic motion 
with increasing pulse length. Two factors contribute to 
this increase. Firstly, a longer pulse obviously means a 
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FIG. 3: Mean displacement of ions (Ar) of neon clusters at 
peak intensity (t = 0) for a full calculation (solid lines) and 
reduced calculations (dashed lines) neglecting field ionization. 
Dependence on the pulse duration T (left, for Neisoo) and the 
cluster size N (right, T = 10 fs), respectively. The average 
(Ar) is shown for all cluster ions in the upper panels and for 
only the inner half of the cluster ions in the lower panels. 



longer propagation of an ion with its acquired momen- 
tum, and secondly the charging of the cluster is higher for 
the longer pulses due to the inherent time scale of about 
2.5 fs set by Auger decay. Although the constant total 
photon number n means that at peak intensity each atom 
will have undergone roughly one photo-ionization event 
for all pulse lengths, clusters under longer pulses will ad- 
ditionally have seen Auger decays. Many of these Auger 
electrons will escape the cluster, cf. Fig.[l] leading to a 
higher charging. The induced ion motion depends also on 
the sample size and shows a strong increase of the mean 
displacement with increasing cluster size (Fig. [3b). Due 
to the additive nature of the Coulomb force a larger clus- 
ter has larger electrostatic energy density than a smaller 
cluster of equal charge density. This leads to stronger 
ionic motion in larger clusters. 

As expected, the displacement is up to 20% lower 
(Fig. [3k) in the full calculation, since the electrostatic 
energy of a charged conductor (with the charge localized 
on the surface) is smaller than that of a homogeneously 
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FIG. 4: Mean displacement (a) and R factor (b) of ions 
of Neisoo cluster embedded in Heisooo droplet (full line) 
and without helium droplet (dashed line) for varying pulse 
lengths. The gray arrow indicates a typical improvement in 
the damage tolerance, see text. 
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charged sphere of equal charge. However, the displace- 
ment averaged over all ion positions does not reflect the 
core-shell nature of the cluster ions with the almost static 
core and the exploding shell in the case of effective field 
ionization. Therefore, Fig. [3]: and Fig.[3ji show the dis- 
placement for the ionic core (the inner half of the ions) 
only, revealing a dramatic effect: The mean displacement 
of the cluster core atoms is reduced by up to 75% for 
Nei5oo for pulses of 20 fs (Fig. [3k) and similarly for the 
larger cluster Nei 50 oo at 10 fs (Fig.[3H) length. Yet, there 
seems to be a lower limit for (Ar) as is most obvious in 
Fig.[3Ji, but also appears through the almost linear in- 
crease with pulse length T in Fig.[3j;. This lower limit 
is due to the recoil acquired by the neon ions during ab- 
sorption of the f2keV photon, which is about O.Of 7 A/fs. 

The different ion dynamics, split into an inert core and 
an exploding shell, naturally leads to the idea of a sacrifi- 
cial layer for imaging experiments as proposed before fToj 
and further investigated in [2l| for carbon-based samples. 
Due to the different material used the results are not di- 
rectly comparable. However, the actual dynamics of an 
exploding shell may be very different in our case since 
field-ionized electrons are created much faster than those 
due to electron-impact ionization from trapped Auger- 
electrons. We consider here a Neisoo cluster embedded in 
a helium droplet of 15000 atoms. During its buildup, the 
net positive charge is efficiently and quickly transferred 
from the cluster to the droplet. The entire neon cluster 
now forms the core of the composite system and remains 
mostly unscathed, while the helium droplet, taking the 
role of the ionized shell, explodes as can be seen in Fig.QJ, 
where (Ar) is shown for the case with and without the 
helium droplet. 

In the context of coherent diffractive imaging it is im- 
portant to ascertain that the strong reduction in (Ar) 
translates into an equally improved quality of the ob- 
tained diffraction pattern. We used our cluster data to 
calculate the diffraction pattern corresponding to a spa- 
tial resolution of 2A, where as in previous investigations 
[Io| we neglected contributions from the tamper and the 
plasma electrons. This diffraction pattern, characterized 



by the intensity /[ eal registered at pixel i of a detector 
having k pixels, is then compared to an ideal diffraction 
pattern /' doal without radiation damage. We define 
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which differs from a previous definition 0] in the nor- 
malization of yj 1} . The R factor of Eq. (JD) therefore 
measures discrepancies due to electron loss (either homo- 
geneous or inhomogeneous) and atomic motion. Both in- 
fluence crucially the image quality which can be obtained 
with XFEL pulses. In comparison to the pure cluster we 
find a substantial reduction of the R factor for the he- 
lium embedded cluster (Fig.|4jD). Much longer pulses are 
tolerated by the embedded system until a similar level of 
damage is reached as in the pure cluster. For example, 
an increase in pulse length by a factor of more than five 
(from T = 2.5 fs to T = 14 fs, cf. gray arrow in Fig.UJa) 
becomes possible by embedding the cluster in helium. 

We have presented detailed investigations for smaller 
clusters which have revealed the importance of internal 
field ionization for subsequent charge migration. For 
larger clusters (R > 60 A) trapping of photo-electrons be- 
comes possible and for clusters with radii in the hundreds 
of Angstrom range, most photo-electrons will be trapped. 
Photo-electron trapping limits the average charge per 
atom in the cluster and the general rule of the surface 
field-strength linearly increasing with cluster size is no 
longer valid. A simple estimate using the Bethe-rule 
for this scenario nonetheless still predicts appreciable 
ionization through quasi-static internal fields. For an 
R = 250 A cluster, for example, in which 50% of the 
atoms are singly photo-ionized, despite the trapping of 
most photo-electrons, the surface-field will still be suffi- 
cient for double ionization of a neutral atom. Hence, the 
field-ionization induced, ultrafast charge migration can 
be expected to play an important role for the success of 
single-shot coherent diffraction imaging experiments for 
a wide array of the samples of interest. 
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